Load data

Summary Heatmaps

By Metabolites

set.seed(194)
suppressMessages(
  dh_heatmap(met_cells, col_names = c('A549_WT', 'A549_KO', 'BEAS2B_WT', 'BEAS2B_OE')))

By Sub-pathway

set.seed(194)
suppressMessages(
  dh_heatmap(met_pth, col_names = c('A549_WT', 'A549_KO', 'BEAS2B_WT', 'BEAS2B_OE')))

A549 xCT KO vs A549 WT

Fold-change plot by metabolites

for (i in unique(met_cells$Super_Pathway)){
  n = 30
  if(length(which(met_cells$Super_Pathway == i)) < n ) {
    n = length(which(met_cells$Super_Pathway == i))
  }
  L2FC_barplot(met_cells, 
               max_items = n, 
               superpathway = i, 
               state_col = 'A549_state', 
               L2R_col = 'A549_L2R',
               sum_subpathway = F)
}

Fold-change plot by sub-pathway

for (i in unique(met_pth$Super_Pathway)){
  n = 30
  if(length(which(met_pth$Super_Pathway == i)) < n ) {
    n = length(which(met_pth$Super_Pathway == i))
  }
  L2FC_barplot(met_pth, 
               max_items = n, 
               superpathway = i, 
               state_col = 'A549_state', 
               L2R_col = 'A549_L2R',
               sum_subpathway = T)
}

BEAS2B xCT OE vs BEAS2B WT

Fold-change plot by metabolites

for (i in unique(met_cells$Super_Pathway)){
  n = 30
  if(length(which(met_cells$Super_Pathway == i)) < n ) {
    n = length(which(met_cells$Super_Pathway == i))
  }
  L2FC_barplot(met_cells, 
               max_items = n, 
               superpathway = i, 
               state_col = 'BEAS2B_state', 
               L2R_col = 'BEAS2B_L2R',
               sum_subpathway = F)
}

Fold-change plot by sub-pathway

for (i in unique(met_pth$Super_Pathway)){
  n = 30
  if(length(which(met_pth$Super_Pathway == i)) < n ) {
    n = length(which(met_pth$Super_Pathway == i))
  }
  L2FC_barplot(met_pth, 
               max_items = n, 
               superpathway = i, 
               state_col = 'BEAS2B_state', 
               L2R_col = 'BEAS2B_L2R',
               sum_subpathway = T)
}